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Abstract 

, We find a family of exact solutions to the semi-classical equations (includ- 

ed ! 

, ing back-reaction) of two-dimensional dilaton gravity, describing infalling null 

OO ; 

matter that becomes outgoing and returns to infinity without forming a black 

as : 

cr 

u '. 



hole. When a black hole almost forms, the radiation reaching infinity in 
advance of the original outgoing null matter has the properties of Hawking 
radiation. The radiation reaching infinity after the null matter consists of a 



brief burst of negative energy that preserves unitarity and transfers informa- 

b : 

tion faster than the theoretical bound for positive energy. 
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Recently we presented a modified two-dimensional (2D) dilaton gravity theory that is 
exactly solvable semiclassically |IJ. In this 2D theory (as in 4D Einstein gravity), infalling 
null matter forms a black hole only if its energy M is above a certain critical value M cr . 
In Ref. U we studied the supercritical case, M > M cr , in which a black hole is formed 
and evaporates by emitting Hawking radiation. In this work we study the subcritical case, 
M < M cr , in which the infalling matter becomes outgoing and escapes to infinity without 
forming a black hole. Because we can obtain the numerical solution to the future of the 
classical outgoing matter, we can see how correlations among outgoing created particles 
preserve unitarity and how Hawking radiation originates as the energy of the infalling matter 
approaches M cr . 

The effective action of the modified 2D dilaton gravity theory is [jl| 
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S = 27 J d2x V~ 9 ^ ~ K ^ R ^ W 

1 N 

+ (4e" 2 * + k)(V0) 2 + 4AV 2 ^ (V/,) 5 

1 i=i 



— (d 2 x^-g{x) [d 2 x'^-g(x')R(x)G(x, x')R(x'), 



8ir 

where is the dilaton field, R(x) is the 2D Ricci scalar, A is a constant, are N matter 
(massless scalar) fields, n = HN/12, and G(x, x') is an appropriate Green function for V 2 (for 

more details see Ref. 0). In null coordinates, z 1 * 1 , and a conformal gauge, g ++ = g = 0, 

g + _ = |e 2p , the equations of motion (from varying and <?+_) take the form d + d_Y = 
= d + d_fi and d + d_X = —X 2 e~ 2Y , where X = e -2 * and Y = — p. The constraints (from 
varying g ++ and g^) are -d 2 ± X - 2d±Xd±Y - (Tl±) c£ + k [(d±Y) 2 + d±Y + t±(^ ± )] = 0, 
where (T±±)c£ = \ Y^i=i{d±fi) 2 is the classical (zero order in h) contribution to the energy- 
momentum tensor of the /j matter fields, and t±(z ± ) are integration functions determined 
by the specific quantum state of the matter scalar fields [Q]. Because is a scalar, but p is 

a function of g^ and d + d-Y = 0, it is possible to define coordinates x + (z + ) and x~(z~) in 

which Y = — p is 0. In this work we use these "Kruskal" coordinates x + ,x~ . 

A general solution in the Kruskal gauge for the X-field equation of motion and constraints 
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IS 



/x + rx^ r f 

d%2 I dxi (Tl + ) c e — Kt+(xi 

— j dx 2 / dx± (t£_)c£ — Kt-(x^) . (2) 



The vacuum solutions are the ones for which (T±±) C £ = and t±(x ± ) = (2x ± )~ 2 They are 
static solutions because in the asymptotically flat coordinates a ± = tier = ±A -1 log(±Aa; ), 
they take the t-independent form X vac (a) = e 2Xcr — nXa/2 + C, where C is a constant. For 
C > C* = (ac/4) (log(«/4) — 1) the vacuum solutions have a null singularity at x^ = 0, 
which is a finite affine distance from any point in the interior of the space-time. For C < C* 
the vacuum solutions have a time-like (naked) singularity on the curve a = a s , for which 
X vac {a s ) = 0. The solution with C = C* is a semi-infinite throat which is everywhere regular 
and geodesically complete ||. 

In addition to initial conditions giving (T+ + )^ and t + {x + ) on asymptotic past null infinity, 
we must also determine (Tl_) c e and t-(x~). To do so, we choose reflecting boundary 
conditions on the matter fields j\ We consider very localized infalling matter for 

which (T+ + ) c ^(x + ) is non-zero only for \x + — Xq\ < e, with Ae << 1 (see Fig. 0). We take 
t + (x + ) = (2x + )~ 2 everywhere, which corresponds to no quantum radiation on Q' - . 

The reflecting boundary conditions are imposed at the boundary curve on which 
X(x + ,x~) = Xb = constant. We take Xb > so that the dilaton field 4> = — |log(X) is 
real. For a dynamical boundary curve, following a general trajectory ) = p{x ), 

the reflecting boundary condition is 



tL = ( P 'M + + K(p'y /2 di( P r 1/2 , (3) 

where ' = d/dx~, and T* is the energy-momentum tensor of the matter fields, including 
both the classical and one-loop contributions. Namely, T' = (T^ l ,) c i + (T^ u ), where (7±±) = 
K>[d±p — (d±p) 2 — t±(a; ± )] is the one-loop contribution to the energy-momentum tensor of 
the scalar matter fields || . The last term on the right-hand-side of Eq. arises because of 
quantum particle creation from the boundary (which is effectively a "moving mirror" 0). 



The classical part of obeys (Tl_)d = (p') 2 (T+ + )rf, which is equivalent to the fields ft 
satisfying either Neumann or Dirichlet boundary condition. While the quantum part obeys 

(T ) = (p') 2 (T ++ ) + Kip'Y^d^L (p')~ l l 2 . The boundary condition (|j) is also conformally 

invariant |7|||. 

To the past of the infalling null matter, x + < Xq — e, i.e., region I in Fig. [I], the geometry 
is one of the static vacuum solutions. The reflecting boundary condition (|3]) implies that 
(T[_) C £(x~) = and t_(x~) = (2x~) -2 for x~ < Xb(x~q — e). The solution (0) to the past of 
the outgoing classical null matter, regions I and II in Fig. [I], is therefore 

X<(x+, x~) = -x + (X 2 x- + P+(x + )) - ^log(-A 2 x+x-) + + C, (4) 
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where M(x + ) = A J x Xi (T+ + ) ct (xi)dxi and P + (x + ) = J x cl {x~i)dx~i are the mass 

and momentum of the classical infalling matter. 

For initial static geometries with C < C*, the (very localized) infalling matter forms 
a black hole only if M > M cr = -X 3 Xq Xq - «A/4 , where M = M(x + > x^ + e) is the 
total mass of the infalling matter, and Xq = x^(xq). The black hole evaporates by emitting 
Hawking radiation and the semiclassical evolution seems to be non- unitary |]||. In 
the subcritical case, M < M cr , the classical infalling matter is reflected from the boundary, 
which is always time-like (see Fig. [l|) and the evolution is unitary. We next find the solution 
numerically to the future of the outgoing classical null matter (i.e., in region III of Fig.|I]) in 
this subcritical case. 

In Kruskal coordinates we use the constraint equations to write the boundary condition 
(0) (on the one loop term) in the form 

/ A1/2Q2 ( v ')-l/2 

*_CO = (pO'VCpGO) - ^- ,; P 9 l , ■ (5) 

1 + k/[2X b ) 

Since we are considering localized infalling matter, the term (p') 1 ^ 2 d 2 _(p')^ 1 ^ 2 , which is of 
the order of M/(Ae 2 ) in the region of classical reflection, is large. However, in this work we 
consider the case in which Xb is small enough |J that XbM/ (/tA 3 e 2 ) << 1. This is consistent 
with the condition, Ti/Xb « 1, that is necessary for the semiclassical approximation to 



be valid everywhere, since in the large N limit we can take % to zero sufficiently fast while 
keeping k = Nh/12 constant [|I]||E§- In the limit X B M/(nX 3 e 2 ) — > 0, the final term in Eq. 
(|) is negligible, and Eq. (Q) reduces to t_(x~) = (p') 2 t + (p(x~)) = {p'/p) 2 /4:. 

One can study the solutions for general infalling matter (not necessarily a very localized 
one). The main features regarding correlations and Hawking radiation with general infalling 
matter |10| are already present in this case with Ae << 1. In this limit the solution to the 



future of Xq (region III in Fig. [T]) can be written in the form 

X>(x + , x~) = -x + (X 2 X - + P+) - ^log(Aa; + ) + F(x~), (6) 

where P + = P + (x + > x^+e) is the total momentum of the classical infalling matter and the 
function F(x~) is determined by the boundary conditions. From and the constraints, 
we see that /tt_(x _ ) = d 2 _F(x~). We define the dimensionless coordinate y = \x~ + P+/X 
and the dimensionless function q(y) = Xxg(x~(y)), and using (|B|) get an ordinary (second 
order) differential equation for the boundary curve Xx~^ = q(y) 

net 

(7) 



Kj d 2 q 

y q{y) + ^ 



dy 2 q{y) 



/ \ 2 

k 1 I dq\ ^dq 



2q 2 (y) \dy) dy 

Define yo = Xxq + P+/X to be the value of y just after the classical reflection (in 
the limit Ae — > 0). Using the continuity of the solution and the classical reflecting 
boundary condition, we get the initial data for q{y): q(yo) = Xxq and (dq/dy)(yo) = 
Xxq M/(Xyo) — Xxq — (4Ax(7)~ 1 /(Xx^yo + k/4). We choose Xx$ = —Xx^ = 1 and 
k = 10~ 6 , and solve (^) numerically for different values of the infalling mass. The re- 
sults that we get are qualitatively the same for any value of k as long as kX << M cr . We 



use an embedded fifth order Runge-Kutta ODE integration routine ||11|| . The solutions for 
q(y) are shown in Fig. |^. 

When the infalling mass, M, is much less than M cr , as in Fig. [2]a, the boundary curve 
is hardly affected by the infalling matter, but when M approaches M cr , as in Fig. ^b, the 
timelike boundary curve is strongly affected and approaches a null curve. 

To see that the solutions are indeed stable (i.e., that the total amount of energy radiated 
to 3 + is finite) we calculate the created energy radiated to S + , (T )(w) = (kA 2 /4)[1 — 



Ay 2 (u)t-(y(u))], where u = — A _1 log(— y) and v = A _1 log(Ax + ) are null coordinate in which 
the metric is manifestly asymptotically flat (i.e., p —>■ on The results are shown in 

Fig. [3] for the same cases considered in Fig. |2|. 

The negative-energy radiation to the future of the original outgoing null matter (i.e., in 
the region u > u ) approaches zero exponentially fast. The magnitude of the total radiated 

negative energy is bounded by | J^(T {u))du\ < (/tA/4)log(4M cr /A/t). If At is the time as 

measured by an asymptotic observer in which this energy E is radiated to S + , then we obtain 
a quantum inequality, \E\At ~ n(M cr -M) 2 (log[M cr /(M cr - M)}f /(MM cr ) < k = Nh/12, 
of the type discussed by Ford and Roman fl2f . The Heisenberg time-energy uncertainty 
principle implies that an attempt to measure the energy of this burst in the available time 
At will disturb the energy by at least \E\. 

We have verified numerically that the magnitude of negative energy radiated to the future 
of the classical reflected matter (u > uq) is equal to that of the positive quantum radiation 
reaching 3 + before the classical reflected matter (u < Uq). Thus energy is conserved. The 
total amount of energy on S + is just that of the original reflected null matter, Eq+ = 
JZo( T --)( u )du = JZo( T --)cedu = M, and the solutions are indeed stable. 

In Fig. [3]a the mass is far below the critical mass and the quantum flux of radiation is 
much less than that of thermal Hawking radiation, while in Fig. |||b the flux of radiation 
before the classical reflection (for u in the range u — A -1 < u < u ) approaches that of 
Hawking radiation. As M approaches the critical mass, the radiation before the classical 
reflection becomes indistinguishable from Hawking radiation originating from a black hole. 

Since the semiclassical evolution of this reflecting solution is unitary by construction 
^3| , [14] j5|j3[] , the radiation to the future of uq must be strongly correlated with the radiation 
to the past of Uq. To see this explicitly we calculate the correlation function 

(x,x') = {T^T^ix')) - (T^(x))(T^ u ,(x')). (8) 

We are especially interested in C : (u, u') which describes the correlations in the outgoing 

radiation on S + . For reflecting boundary conditions there is a closed form expression 0] : 
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1 {d u v B {u)f (d u ,v B {u' )f 
Stt* Mu )-v B (u')) 4 



— K ) = ^ — ; — — r^Ti — > ( 9 ) 



where v B (u) is the boundary curve in the coordinates v = A _1 log(Ax + ) and u = 
— A _1 log(— (\x~ + P + /X)). Let us first consider the correlations in the radiation reaching 
$5 + before the classical reflected matter. For u < uo we have v B (u) = A -1 log(— X 2 Xq Xq ) — 
A -1 log(e _A " + P + /X), and for M near M cr we have e Xu ° » 1. Taking u < Uq and vl < Uq, 



with both near uq, we find to leading order in Ae Xu /P + and in Ae Au '/P_ 



+ • 



A 4 fe- Au ) 2 fe- Au,N ) 2 / Ae" Au \ 
C (u, u') = Ar f t M \ L 1 + 0(— ) . (10) 



The leading term in (|T0f) is just the correlation function in 2D for thermal black body 
radiation 0] with temperature T = X/2tt, the temperature of 2D dilatonic black holes 
T5| , [T6|| . Through the same calculation that gave Eq. (|ToD, we find that fllCf ) is also the 



correlation function for Hawking radiation from an evaporating 2D black hole. 

To explicitly verify the correlations between the quantum radiation reaching 9 + before 
and after the classical reflected matter, we numerically calculate @ by holding u' fixed on 

either side of uq and plotting C i (u, it') as a function of u. The results are shown in Fig. 

^, where the mass of the infalling matter is just below the critical mass, M = 0.99M cr . 

In Fig. ^a we show the correlations between the radiation at Xu' = 3.6, and the radiation 
elsewhere on £y + . This value of Am' is sufficiently close to Xuq(= 4.6) that the radiation has 
the same form as the Hawking radiation from a black hole. For u < uq the correlation is 
almost thermal, as in Eq. ([10|) . It diverges at u — u' and becomes small for \u — u'\ ~ 1/A. 
The most dramatic effect occurs just after the reflected matter reaches S + , i.e., at u = u + 5. 
The correlations rise (continuously when e is finite fllQ| ) to extremely high values of order 
e 15 . These huge correlations are required for unitary evolution because the negative-energy 
radiation reaching 3 + after the classical reflected matter is very localized in time. 

Fig. |]b shows the correlations calculated numerically between the negative-energy radi- 
ation reaching S + just after the classical reflected matter and the radiation elsewhere on 
In this figure we take Am' = Xuq + 5 and 5 — > + . The correlations diverge at u = u' 



and decrease rapidly when u > u' . On the other hand the correlations with the Hawking 
radiation (u < Uq) remain extremely high for relatively large values of X\u — u'\. In contrast, 
if one were to take An' >> Xuq then the only significant correlations would occur for Xu 
near Xu'. 

The above results may plausibly be interpreted as arising from the creation of particle- 
antiparticle pairs []17| |. The particles reach infinity and give rise to the positive (Hawking) 



radiation, while the antiparticles carrying negative energy are reflected from the boundary 
and give rise to the negative energy radiated to S + after the classical reflected matter. If it 
were not for the negative-energy burst the correlations between the particles and antiparticles 
would be lost and the final state would be a mixed state. Thus, the negative-energy burst can 
be regarded as carrying information equal to the magnitude of the entropy S of the Hawking 
radiation reaching S + ahead of the classical reflected matter. As M approaches M cr , where 
M cr » k\/A, we find that S « (iV/12)log[4M cr /(«A)]. The negative-energy burst thus 
carries information to S + at the rate / = S/At. Using the earlier upper bound on \E\At, 
we find that this gives / > (iVA/48)(i — M/M cr )~ 2 . The theoretical upper bound on the 
bulk rate of linear information flow in iV channels for E > is (TTSf I max < [E/ (27rh)]\og 2 {N) 



to within a factor of order 1. If this bound can be extended to E < by replacing E by \E\, 
then it gives for the present system: I max < [A7"A/(967r)]log[4M cr /(reA)]log 2 (iV). This bound 
is clearly exceeded for large but finite N in our system as M approaches M cr . 
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supported by the National Science Foundation under grant PHY 95-07740. 
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FIGURES 

FIG. 1. Penrose diagram for a typical subcritical solution. The 2D space-time is the inte- 
rior of the ID null hypersurfaces 9 + and the time-like boundary curve, which intersect at 
%~ = (x + = 0,x~ = — oo), i° = (x + = co,x~ = — oo) and i + = {x + = oo, x~ = 0). 

FIG. 2. The boundary curve Ax^ = q(y). In (a) the mass of the infalling matter is M = 0.1M cr 
and in (b) it is M = 0.9M cr . The dimensionless coordinate y = \x~ + P + /\ is plotted on the 
horizontal axes. 

FIG. 3. Quantum radiation on 9f+. In (a) M = 0.1M cr and in (b) M = 0.9M cr . Here (T__) is 
in units of kA 2 /4 (which is the value for Hawking radiation from 2D black holes). In (a) Auo = 0.1 
and in (b) Auo = 2.3. 

FIG. 4. Log of the correlation function C , (u,u r ) for fixed u'. In (a) we take An' = 3.6, 

and in (b) Xuq approaches Xuq = 4.6 from above. 
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FIGURE 1 
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